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We study the magnetic properties of two types of one dimensional XX spin 1/2 chains. The first 
type has only nearest neighbor interactions which can be either antiferromagnetic or ferromagnetic 
and the second type which has both nearest neighbor and next nearest neighbor interactions, but only 
antiferromagnetic in character. We study these systems in the presence of low transverse magnetic 
fields both analytically and numerically. Comparison of results show a close relation between the 
two systems, which is in agreement with results previously found in Heisenberg chains by means of 
a numerical real space renormalization group procedure. 
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One-dimensional quantum spin systems have been ex- 
tensively studied over the last yearsi. In particular, ran- 
domness has a profound effect on their physical proper- 
ties and is always present in real systems through impu- 
rities or structural disorder. This can even produce sin- 
gular behaviors in the magnetic properties, not observed 
in pure systems. One of the main motivations to study 
disordered quantum spin chains is the possibility of clas- 
sifying their behavior in universality classes associated to 
different regions in their phase diagrams**. 

In the last few years, numerical works^* 3 - have shown 
that the thermodynamic properties of disordered Heisen- 
berg chains with nearest neighbors (NN) and next near- 
est neighbors (NNN) couplings, both antiferromagnetic, 
are very similar to those found in disordered chains with 
only NN couplings which can be either antiferromagnetic 
or ferromagnetic. In fact, it was shown that under Real 
Space Renormalization Group (RSRG) the former sys- 
tems flow to a fixed point characterized as a chain with 
only NN couplings in a given distribution, taking both 
antiferromagnetic and ferromagnetic values. 

More specifically, let us consider a Heisenberg chain 
with Hamiltonian 



H 



i+2 



(1) 



where Si are spin-1/2 operators and the couplings 3{ > 
and J[ > 0, both antiferromagnetic, follow probability 
distributions P(Ji), P{J[)- Let us review the arguments 
of Ref. • K we consider the adjacent spins that are 
coupled by the strongest bond (say the spins 3 and 4 in 
the figured]), and its neighbors, we have a problem whose 
Hamiltonian can be written as 



H = H + Hj + H rest , 



(2) 
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FIG. 1: Schematic picture of a disordered-antiferromagnetic 
NN-NNN chain. 



where 



H 



J34S3 ■ S4 



Hi — J23S2 ■ S3 + J45S4 
+ J24S2 • S4. 



S, 



J35S3 ■ £5 



(3) 



(4) 



and H rest corresponds to all the other spins that are not 
coupled to spins 3 and 4. 

The sole consideration of H + Hi is enough to deter- 
mine the effective interaction between the spins 2 and 
5, as follows^^. The ground state for the Hamiltonian 
Hq is a singlet with energy Eq = —ft/34 while excited 
triplet states have energy E\ = Since J34 is the 

largest bond, one can take Hi as a perturbation. Re- 
garding S2 and S5 as external operators, a second-order 
perturbation calculation gives an effective Hamiltonian 
describing the low energy sector (we consider only the 
coupling generated between 2 and 5) 



E = 
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From this result one can remove the spins S3, S4 from 
the original Hamiltonian, replacing them by an effective 
NN coupling J25 



■I 



25 



(J23 " J24XJ45 - J35) 

2 J34 



(6) 



uniform magnetic field. By means of the Jordan- Wigncr 
transformation^ 



(8) 



between S2 and S5. 

Following this RSRG decimation procedure one ends 
up with a NN spin 1/2 chain in which effective interac- 
tions like J25 can be ferromagnetic. Notice that when 
the NNN couplings are very weak compared to the NN 
couplings the ferromagnetic effective interaction is un- 
likely to appear. On the contrary, for strong NNN cou- 
plings, RSRG does generate effective ferromagnetic cou- 
plings and the system flows to a phase controlled by large 
effective spins at low energies^. 

It is the aim of the present note to further understand 
the connection between these and related systems, in 
more general situations. In particular we are interested 
in comparing magnetic properties of easy plane XX spin 
1/2 chains in the presence of uniform magnetic fields. 

We first study analytically the magnetic properties 
of quantum XX spin 1/2 chains with antiferromag- 
netic and ferromagnetic NN interactions, using an el- 
egant argument relating magnetization with random 
walk problems^. In a second step we implement a nu- 
merical self consistent method based in a mean field 
approximation^, which allows for analyzing quantum 
XX spin 1/2 chains with NN and NNN interactions, to 
be applied in the antifcrromagnetic case. Finally, com- 
mon features found in these spin systems are discussed. 

The structure of the paper is the following. Analyti- 
cal results derived from the random walk argument arc 
presented in section [n] for homogeneous and dimerized 
disorder in antiferromagnetic NN chains, as well as novel 
results for dimerized distributions of antiferromagnetic 
and ferromagnetic couplings in NN chains. The self con- 
sistent numerical method is presented in Sect ion [Hi] then 
tested on pure (ordered) NNN antiferromagnetic chains 
in Section HVl and finally applied to the NNN disordered 
case in Section [V] Numerical results are compared with 
exact diagonalization for small chains, up to 24 spins, at 
each step. The comparison between both types of chains, 
summary and conclusions are presented in Section [VII 



II. ANALYTICAL TREATMENT FOR NN SPIN 
1/2 CHAINS 



In this Section we study analytically the magnetization 
of one dimensional spin 1 /2 systems with NN interactions 
whose Hamiltonian in the XX model is 

H = E Ji(SfSf +1 + SfSf +1 ) -h^Sf, (7) 



where Jj are random nearest neighbor couplings (cither 
antiferromagnetic (AF) or ferromagnetic (F)) and h is a 



where 4>i = X);=i c l c i> we can rewrite this Hamiltonian 
in terms of spinless fermionic operators 



H = J2 U (4^+1 + H.c.) - 



(9) 



(in this Section we use t = J/2). Notice that the mag- 
netic field acts as a chemical potential for the fermions. 

We start with the study of the homogeneous disordered 
case where the couplings follow a homogeneous distribu- 
tion {P(ti) = P{tj) Vi,j). Following the argument pre- 
sented in Ref. @ one can introduce a random variable Ui 
which, for a one particle cigenstatc of the Hamiltonian 
in eq. (|9|), undergoes a random walk behavior between 
a reflecting barrier at u max = \og{t 2 / E) and an absorb- 
ing barrier at u m i n = log(E) with E being the energy 
of the eigenstate and t the positive average value of AF 
couplings ti. The relevant quantity to compute here is 
the number of eigenstates N(E) with energies below E, 
which for low energies is related to the average number 
n of steps necessary to complete a diffusion cycle from 
the reflecting barrier to the absorbing one. The relation 
between them is^ 



N(E) = — + -. 
y ' 2n 2 



(10) 



In the present case one gets n ~ {u max — u m i n ) 2 /a 2 , 
where a 2 is the variance of the coupling distribution P, 
rendering 



N(E) 



1 



1 



{\n{t/E) 2 f 



(11) 



This result indeed describes the magnetization of the spin 
system ([7]) . We write the mean magnetization m in terms 
of fermionic variables using the Jordan- Wigncr transfor- 
mation 



<E^HE« c U)-l/2), 



(12) 



in order to exhibit the relation between magnetization 
and fermion filling. Recalling that the magnetic field 
acts as a chemical potential and regulates the fermion 
filling, for low magnetic field we finally obtain a mean 
magnetization per site 



M(h) ~ 1 ( 5 2 
{> 2 \[hx(t/h) 2 ] 2 



(13) 



(here M = 2m/N, with N the number of spins in the 
chain, is normalized to 1 at saturation). 
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Now we turn to the case of dimcrized inhomogencous 
distributions ((P dd(ti) f° r °dd sites i different from 
Peven{tj) for even sites j). The result above was suitably 
generalized in Ref. for this case: the random variable 
Ui undergoes a random walk with diffusion coefficient D 
and a driving force F given by 

D = i[var2 rfd (log(^)) + 2var^ en (log(i J 2 ))], (14) 
F = (log(t 2 )) odd - (log(if )) ewen , (15) 

where var stands for the variance of the correspond- 
ing distribution. The average numbers of steps for a 
diffusion cycle to be completed behaves now as ri ~ 
e" ( u "x" -tl ™in)/ 2 . Then the magnetization of the system 
is seen to follow a power law 

M(h) ~ h a , (16) 

with a = 

In what follows we generalize this procedure to study 
the system of interest here, namely a disordered spin 1/2 
chain with AF and F NN interactions. Let us consider 
the following binary coupling distribution 

Pfa) = xP F (ti) + (1 - x)P AF (ti) > (17) 

with weight x for F couplings and 1 — x for AF ones, 
combined with dimerization in the sense described above 
(both Pp and Paf ar< 3 different for even and odd sites). 
A similar pattern for disorder was proposed in recent 
numerical studies 2 ^. 

We can again map the system onto a random walk 
problem with driving force and appropriate barriers. In 
particular the driving force can be written in terms of 
the single distribution parameters as 

F = x(log((t iodd ) 2 /(t ieven ) 2 )) F + 
(l-x)(log((t iodd ) 2 (t ieve J 2 )) AF , (18) 

where subindcxes odd and even indicate the distribution 
to be used for disorder average. 

Notice that even under the strong hypothesis that both 
single distributions are dimerized (inhomogeneous) 

(log((( !o J 2 /(«w) 2 )>^0, 

and 

(log((t iodd ) 2 (t ieven ) 2 )) AF ^0, 

the competition between AF and F couplings can even- 
tually cancel the driving force, under the condition 

{\og{tlJtl_)) AF = (19) 
x[(\og(tlJtl v J) AF - (\og(tl d Jtl v J) F }. 

This shows that there are two phases present in the sys- 
tem. For the coupling distribution in eq. (jTTJ) we have 



that, at least for low magnetic field h, the magnetiza- 
tion follows a power law M ~ h a in most of the pa- 
rameter space, while there exists a line, with x satisfying 
cq. (|20p . where the magnetization is logarithmic in h, 

M ~ iog(?i^)^ ■ 

It is important to stress that in the power law regime 
the dynamical exponent a can be larger or smaller than 
one. For the present case it is still given by 2F/D where 
F and D are computed as in eqs. (|15p and ([M]) but with 
the binary even and odd distributions in eq. (|17p . When 
the disorder parameter (in this case var (ln(t 2 ))) is small 
and the dimerization is no longer considered, the dynami- 
cal exponent turns out to be larger than one. In contrast, 
when the disorder is stronger the dynamical exponents 
take values smaller than one and the magnetic suscepti- 
bility displays a singularity at the origin. It is worth to 
point out that this behavior was also reported by RSRG 
studies in dimcrized NN chains^. 



III. NUMERICAL SELF CONSISTENT 
TREATMENT FOR ANTIFERROMAGNETIC 
NNN SPIN 1/2 CHAINS 

In this Section we consider a spin 1 /2 XX Hamiltonian 
in d = 1 with both NN and NNN AF position dependent 
interactions under a uniform magnetic field, namely 

N 

H = ^(-Sf^i + W+i) 

i=l 

N 

+ J[(S*S* +2 + SfSf +2 )} - h J2 St , (20) 

i=l 

where N is the number of spins in the chain and periodic 
boundary conditions are assumed. The inclusion of NNN 
couplings does not allow the analytical procedure used 
above. We then perform a numerical self consistent mean 
field (SCMF) study of this system. 

We first review the procedure proposed ir*£ and then 
apply it to our present case. In terms of the fermion 
operators in eq. ([5]) this Hamiltonian reads 

N J 

+ Ey[ e_M " I+lc ^ + 2+H.c.]+ (21) 

N 1 

i=l 

First, in order to enable a single particle treatment of H, 
we approximate the local fcrmionic occupation numbers 
n, = c\ci by their expectation values in an arbitrarily 
chosen initial state to be varied and determined self con- 
sistently The local parameters satisfy the constraint 
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~^2iLi( n i ~~ 1/2) = m , with m the system magnetization. 
Then, in this mean field (MF) approximation the Hamil- 
tonian can be written as a quadratic form 



(22) 



where 



Jij({ni}) 



Ji 
2 

Ji 



if i,j are NN, 



i e «tn+i if are NNN, 



, otherwise. 



We have omitted the Zeeman term h^ i=1 {c\ci — i) as, 
being diagonal, it can be added later on. 

It is our aim to find an approximation to the actual 
ground state (GS) at a given magnetization m. Thus, 
we estimate this state in the MF Hamiltonian ([22jl by 
solving the one particle spectrum and filling the lowest 

energy-levels to satisfy the constraint X^ti( n « — V^) = 
to. Then, we compute a new set of local parameters n[ = 
(GS\c\ci\GS) which we use again as input for Eq. (|2"2"j) . 
Iterating this procedure, we finally obtain a fixed point 
configuration of occupation numbers ^({rip}) = n p . 

Specifically, the quadratic Hamiltonian can be written 
in diagonal form 



N 



H = Y^<k)d{d k , 



(23) 



where the operators Cj are related with dk by the unitary 
transformation 



(24) 



where Qik is the matrix of eigenvectors of Jij({ni}). 
Using standard methods^ we can easily compute the 
eigenvalues e(k) of H^lf^ ({rij}) and eigenvectors Qik for 
fairly large spin systems. The set of dk satisfy fermion 
anticommutation relations {dk, dt,} = 5k,k' and the total 
fermion number is conserved 



TV N 

i=i fe=i 



(25) 



So, the ground state with magnetization to, in the diag- 
onal basis is given by 



m+JV/2 

\gs)= n 4io). 
fc=i 



(26) 



In this state the energy at zero temperature (for zero 
magnetic field) is simply given by 




h/J 

FIG. 2: Magnetization curves for pure chains with L = 102 
and L = 84 sites both with J' /J = 0.6. For M = there is a 
magnetic plateau. 



E GS (m,0) 



m+N/2 

E 

fc=i 



e(k). 



(27) 



whereas for h ^ the total energy is shifted as 
Ecsim, h) = Ecsijn, 0) — hm. The magnetization curve 
m(h) can finally be obtained by minimizing the energy 
£"g5(to, h) for different given magnetizations. 



IV. RESULTS FOR THE PURE CASE 

We have tested the SCMF procedure in a pure (or- 
dered) J — J' XX model using chain lengths of up to the 
order of 100 spins, and compared its results with those 
found by exact diagonalization in smaller chains. In Fig. 
[21 we show the magnetization behavior under a magnetic 
field for J' I J = 0.6 which clearly indicates the existence 
of a magnetization plateau at M = 0. 

This initial plateau shows up only in a narrow region of 
J' I J which in the SCMF approximation was estimated 
within the bounds 0.55 < J'/ J < 0.75. No subsequent 
plateaus were observed in the system. Notice that in 
other models (e.g the XX Z modet^) there is a plateau at 
M = 1/3, however this not the case in the XX situation. 

To lend further support to our SCFM approach, 
we compared the above results with those obtained in 
smaller systems using exact diagonalization^. In Fig. [3] 
we exhibit the magnetization curve obtained by Lanczos 
technique using chains of 12, 18 and 24 spins, where an 
M = plateau also emerges at J' /J = 0.6. The regime 
where this plateau appears turns out to be slightly higher 
than that found with SCMF. This can be observed in Fig. 
[4] where we show the magnetic phase diagram obtained 
in a wide region of coupling parameters. Here, each line 
stands for a critical field above which the magnetization 
is increased by flipping one spin. In agreement with our 
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h/J 

FIG. 5: Magnetization curve for a disordered chain with 
L = 102 sites and J'/J = 0.6, averaged over 150 samples. 
FIG. 3: Magnetization curve for small pure chains with The magnetization drops quickly to zero rendering a singular 
J'/J = 0.6, obtained by exact diagonalization. There is a magnetic susceptibility at zero field, as shown in the inset, 
clear plateau at M = 0. 
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FIG. 4: Magnetic phase diagram for pure chains. In the 
regime 0.5 < J'/J < 0.75 the critical field h c increases, show- 
ing the formation of the M = plateau. 

SCMF expectations, we can see that the first critical field 
(lowest line of each studied length) is higher in the region 
0.5 < J' I J < 0.75 where the M — plateau is favored 
and size effects become substantially reduced. 

V. RESULTS FOR THE DISORDERED CASE 

We now apply the SCMF procedure to our main in- 
terest, namely disordered antiferromagnetic XX chains 
with a Hamiltonian given by cq. (|20p . In this case Jj 
and J[ are random NN and NNN coupling exchanges re- 
spectively. For concreteness, let us consider Gaussian 

distributions of exchanges P{Ji) oc e ^ with mean 
value J = J > for NN couplings, J = J' > for NNN 



couplings and the same disorder strength a in both cases 
(we explicitly forbid negative couplings in the Gaussian 
tail). We compute numerically the chain magnetization 
by averaging over many disorder realizations. Typically 
we considered over 200 samples with periodic boundary 
conditions and random sets of initial fermionic distribu- 
tions. 

We focused particular attention on low magnetic fields, 
so as to detect possible singularities near h = as those 
observed in disordered NN chains^. In Fig. [5]wc display 
the magnetization curve for 102 spins with J'/J = 0.6 
and a I J = 0.5. Notice that the magnetization plateau 
observed in the pure case {a = 0) is now totally sup- 
pressed. In fact, all studied values of a suggest that the 
plateau of the pure system is unstable under disorder. 
This observation was also corroborated in smaller sys- 
tems after diagonalizing them exactly over 100 disorder 
realizations. In Fig. [6] we show for comparison the mag- 
netization curves obtained for the same values of J, J' 
and a. On the contrary, the magnetization remains fi- 
nite and drops quickly to zero at zero field. Actually, the 
magnetic susceptibility x = exhibits a divergence at 
h = , as shown in the inset of Fig. [H 

A thorough exploration of mean values of NN and NNN 
couplings and disorder strength shows that the low field 
magnetization curve shows a behavior compatible with a 
power law M(h) ~ h a in most of the parameter space, 
except on a small region J'/J < 10 -4 where M decreases 
in a logarithmic form M ~ (i og {h' 2 ))' 2 ■ ^ n ^ ne power law re- 
gion an exponent a < 1 is obtained for disorder strength 
a I J > 0.3, corresponding to a singularity in the zero field 
magnetic susceptibility. For a/ J < 0.3 the magnetization 
decreases with an exponent a generally larger than one. 

In Fig. [7] we show a schematic diagram of our SCMF 
numerical results. The dashed zone denotes a logarithmic 
behavior, the light gray region represents a power law 
behavior decrease of M with exponent a < 1, whereas in 
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we have found three phases in these systems, the domi- 
nant one being a power law dependence of the magnetiza- 
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FIG. 6: Exact magnetization curves of disordered chains av- 
eraged over 100 samples for L = 24, 18 and 12 spins (solid, 
dashed and dotted lines respectively; notice that averages 
were taken for the critical fields, not for the allowed magneti- 
zation values, so the curve is not smoothed). J, J' and a are 
taken as in Fig. [5] The inset denotes the standard deviations 
a h of the corresponding critical fields. 



FIG. 7: Schematic diagram of the low field results obtained 
for disordered NNN chains. In a tiny (dashed) region the 
magnetization is logarithmic like. In contrast, in most of pa- 
rameter space the magnetization follows a power law, with 
larger disorder leading to singular susceptibility (light gray 
region) and small disorder leading to smooth magnetization 
(gray region). A plateau at M — is present only in the 
absence of disorder. 



the gray zone this exponent is larger than one. For the 
ordered case a = we show a bold line representing the 
magnetic plateau at M = within 0.55 < J'/J < 0.75 . 

We can finally compare the low field behavior of disor- 
dered antiferromagnetic NNN chains with those found in 
disordered NN AF-F chains (Section II). We can stress 
that both systems have two phases: a dominant one char- 
acterized by a power law magnetization behavior and a 
small region of the parameter space where that behavior 
is close to a logarithmic type. Moreover, in the power 
law regime both systems can develop singular or smooth 
zero field magnetic susceptibilities. 

VI. SUMMARY AND CONCLUSIONS 

In the first part of this work we have studied analyti- 
cally XX spin- 1/2 chains with random NN interactions 
both antiferromagnetic and ferromagnetic, by suitably 
extending the analysis based on a random walk prob- 
lem presented in Ref. [||. We have distinguished three 
regions in the parameter space with different low mag- 
netic field behavior: singular power law M cx h a with 
a < 1, smooth power law with a > 1 and a logarithmic 
dependence cx n^n^syp ■ The second part discusses an al- 
ternative approach to random antiferromagnetic spin-1/2 
chains which allows for NN and NNN interactions, using 
a numerical self consistent mean field method adapted 
from Ref. 0. As for the NN chains analyzed before, 



tion with the low external field in most of the parameter 
space. Also, a logarithmic dependence of the magnetiza- 
tion was found just within a small region of parameters. 

Our results show that the two systems have similar 
low energy properties. A similar analogy has been found 
in the study of the same two systems for the Heisen- 
berg (5(7(2)) version^. In both of them there are three 
phases and most of the parameter space is dominated 
by a power law dependence of M{h). Both systems also 
show a small region where the magnetization displays 
the same kind of logarithmic singularity. In both sys- 
tems the power law region follows a dynamical exponent 
a > 1 for weak disorder, thus yielding a well behaved 
magnetic susceptibility. On the contrary, strong disorder 
yields exponents smaller than one and consequently the 
susceptibility becomes singular at h = 0. This later be- 
havior has been also found in dimerized spin-1/2 chains 
with disordered NN frrteractionsii. 
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